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Abstract . 

Subroutine  T"STADETr  analyzes  a  structure  to  find  if 
it  is  stable  and  determinate.   The  structure  must  first  be 
idealized  as  a  rigid  body  supported  by  a  system  of  links. 
The  program  checks  on  the  arrangement  and  number  of  the 
support  links.   Flags  are  set  so  that  this  subroutine  may 
be  used  to  select  whether  to  use  simple  statics  or  strength 
of  materials  methods  in  subsequent  analysis  of  the  problem. 
Testing  of  the  subroutine  is  discussed. 

Two  input-output  programs  are  also  described.   One 
calculates  the  reactions  in  the  equivalent  support  linkages 
from  the  resultant  load;   the  other  assumes  that  the  pro- 
perties of  the  links  are  known  and  calculates  the  allowable 
load. 
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1.   Introduction. 

Digital  computers  have  been  extensively  applied  to 
the  solution  of  mechanics  problems .   However,  these  ap- 
plications have  been  limited  in  scope.   The  computer  has 
been  used  as  a  "super  slide  rule"  to  solve  a  particular 
type  of  problem  which  was  too  time-consuming  or  too  dif- 
ficult previously.   It  is  the  intention  of  Professor  Gilles 
Cantin  to  attack  the  problem  of  the  design  of  structures 
by  a  digital  computer  in  a  general  rather  than  a  specific 
manner. 

One  of  the  first  routines  needed  for  the  Cantin  Pro- 
ject was  a  routine  to  reject  unstable  structures  and, 
perhaps,  make  some  selection  as  to  the  methods  of  analysis 
to  be  used  in  subsequent  tests.,   The  subroutine  T'STADET" 
analyzes  an  idealized  structure  to  find  if  it  is  stable 
and  determinate.   It  sets  flags  which  can  be  used  to  print 
out  a  description  of  the  kind  of  problem  involved.  When 
the  subroutine  finds  the  first  condition  which  would  pre- 
clude the  solution  of  the  problem  by  simple  statics,  a 
flag  is  set  and  control  returns  to  the  input-output  program 
If  the  "error"  is  caused  by  a  redundant  support  system, 
other  methods  of  analysis  may  be  available  whereas  there 
is  no  solution  for  an  "unstable"  problem. 

The  definition  of  stability  which  is  considered  here 
is  more  restrictive  than  the  question  of  stable  equili- 
brium usually  encountered  in  basic  statics  courses.   Here 
is  meant,  not  the  criterion  of  energy  change  with  small 


perturbations,  but  the  geometrical  stability  of  infinitely- 
strong  members.   As  expressed  by  Kinney, 

A  stable  structure  will  support  any  conceivable 
system  of  applied  loads,  resisting  these  loads 
elastically  and  immediately  upon  their  applica- 
tion the  strength  of  all  members  and  the  capacity 
of  all  supports  being  considered  infinite. 

Consider  the  two  links  supporting  a  weight  shown  in 

Fig.  1. 
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Fig.  1.   Support  links  before  and  after  displacement. 

In  Fig.  la,  with  the  links  horizontal,  there  is  no 
vertical  component  of  the  tension  in  the  bars  to  support  W 
Even  with  infinitely  strong  foundations  and  members,  the 
geometrical  arrangement  precludes  a  vertical  component  un- 
til some  downward  displacement  of  W  has  occurred  (Fig.  lb), 
The  methods  of  strength  of  materials  are  required  to  find 
the  equilibrium  position  of  the  structure.   Since  W  is  not 
resisted  "immediately  upon  its  application,"  this  arrange- 
ment fails  the  test  for  stability. 

Both  a  structure  and  its  support  system  must  be 
stable.   The  problem  which  is  examined  here  is  limited 
to  that  of  the  external  stability  of  a  structure.   The 


J.  Sterling  Kinney,  Indeterminate  Structural  Analysis 
p.  20,  1957. 


structure  itself  is  considered  to  be  a  rigid  body  so  that 
only  the  system  of  supporting  links  is  of  concern.   No 
loads  are  imposed  during  the  stability  analysis  since 
equilibrium  is  not  involved . 

However,  in  some  cases  forces  such  as  gravity  must  be 
considered  in  order  to  have  a  structure  which  meets  this 
restricted  definition  of  stability.   This  is  a  familiar 
assumption  in  basic  statics  where,  for  example,  the  gravity 
force  provides  a  reaction  against  upward  motion  of  a  roller 
nest.   It  is  up  to  the  user  of  this  program  to  decide 
whether  a  "gravitational"  restraint  is  sufficient  or  whether 
a  positive  restraint  is  needed. 

The  user  must  first  idealize  his  structure  so  that  the 
subroutine  can  operate  on  it.  Appendix  A,  which  shows  the 
linkage  diagrams  of  some  common  supports,  may  be  helpful  in 
this  regard.  However,  it  is  at  this  point  that  engineering 
judgment  is  required.  Needless  to  say,  if  the  proper  equi- 
valent support  system  is  not  employed,  no  reliance  can  be 
placed  on  any  answers  which  are  produced . 

Fig.  2  shows  the  decisions  which  are  made  by  the  sub- 
routine.  The  titles  within  the  blocks  indicate  the  basic 
sections  of  the  subroutine.  These  titles  are  roughly  indica- 
tive of  the  purpose  of  the  sections  though  more  information 
is  usually  produced  than  the  title  would  lead  one  to  believe. 

For  input  the  subroutine  requires  the  direction  cosines 
of  the  support  linkages,  the  coordinates  of  the  points  of 
application  of  the  links  and  the  number  of  links.   Appendix  H 


Figo  2  Simplified  Flow  Chart 
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shows  how  the  input  must  be  arranged  to  be  read  in  by  the 
input-output  programs  of  Appendix  C.  The  output  consists 
of  the  two  flag  arrays,  KPROB  and  NOSOL. 

The  subroutine  first  considers  the  points  of  appli- 
cation of  the  supporting  links „   If  these  points  are  along 
a  line,  the  question  of  whether  it  is  a  two-  or  three- 
dimensional  problem  depends  on  the  orientation  of  the  links. 
If  the  points  are  in  a  plane,  the  links  must  all  be  in  that 
same  plane  for  a  two-dimensional  problem  to  be  involved.   If 
the  points  are  randomly  located  in  space,  only  a  three- 
dimensional  problem  can  be  involved.   The  two-dimensional 
problem  requires  three  linkages  for  stability;   the  three- 
dimensional  problem,  six.   If  there  are  fewer  linkages,  an 
"unstable"  error  print-out  is  made;   if  more,  an  "indeter- 
minate" one  is  made. 

When  it  has  been  determined  that  the  correct  number  of 
support  linkages  is  involved,  the  orientation  of  the  links 
is  considered .   The  following  restrictions  apply  to  the 
three-dimensional  case: 

1.  All  the  links  may  not  be  parallel. 

2.  They  may  not  lie  in  parallel  planes. 

3 .  The  axes  of  the  links  may  not  intersect  one 
straight  line  since  limited  rotation  would 
be  possible  about  that  line. 

4.  The  axes  of  the  links  may  not  intersect  in  a 
point  or  limited  rotation  is  possible  about 
that  point . 

5.  Four  of  the  six  links  may  not  be  parallel  or 
intersect  in  any  one  point  for,  in  these  cases, 
a  straight  line  can  be  found  through  the  point 
(or  infinity)  which  intersects  the  axes  of  the 

two  remaining  links  and  restriction  3  •>  is  violated. 

In  the  two-dimensional  case  only  restrictions  1.  and 
4°  apply . 


Where  the  links  are  all  parallel  to  a  given  line, 
limited  motion  is  possible  in  a  direction  perpendicular 
to  this  line. 

In  Fig.  3 ,  where  the  pairs  of  links  are  in  parallel 
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Fig.  3«  Arrangement  with  links  in  parallel  planes. 

planes,  some  motion  is  possible  in  a  direction  perpen- 
dicular to  these  planes  before  a  restoring  component  is 
developed.   In  Fig.  4  some  rotation  about  the  intersection 
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Fig.  4.  Arrangement  with  intersecting  lines  of  action 
of  links. 


of  the  lines  of  action  of  the  links  is  necessary  before 
a  restoring  couple  is  developed. 

The  following  section  and  the  appendices  describe 
the  "STADET"  subroutine,  its  auxiliary  subroutines  and 
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the  Input-Output  Programs  which  were  used  to  test  "STADET". 
Fortran  language  was  used  for  all  the  programs .   The  testing 
was  done  on  a  Control  Data  Corporation  1604  digital  computer, 


2.   Method. 

Subtitles  indicate  the  sections  of  the  flow  charts; 
(Figs.  5  through  9)  and  the  "STADET"  list  (Appendix  J) 
under  discussion.   Circled  numbers  on  the  flow  charts  are 
used  for  two  purposes.   The  first  is  to  facilitate  the 
comparison  with  the  list  of  Appendix  J;   the  second  is  to 
indicate  a  transfer  to  another  section  of  the  program 
when  the  use  of  a  line  to  show  the  path  would  be  cumber- 
some . 

The  subroutine  only  operates  up  to  the  detection  of 
the  first  error .  The  appropriate  flag  is  then  set  and  a 
return  is  made  to  the  Input-Output  Program. 

a.   C0PLANAR1 

This  section  of  the  subroutine  generates  the 
position  vectors  which  give  the  locations  of  the  points 
of  application  of  the  links.   It  checks  to  see  whether 
these  points  are  coincident,  lie  on  a  line,  are  in  a  plane 
or  are  randomly  located  in  space.   In  the  case  where  the 
points  of  application  are  in  a  plane,  this  section  of  the 
program  checks  to  see  if  the  links  are  in  the  same  plane 
and  so  determines  whether  a  planar  problem  is  involved. 

After  deciding  whether  a  two-  or  three-dimensional 
problem  is  involved,  this  section  counts  the  number  of  links 
to  find  out  if  the  number  is  sufficient  for  stability. 

The  computer  operates  with  a  binary  approximation  to 
a  decimal  number  when  it  is  used  in  the  floating  point 
mode.   Each  calculation  is  approximate  and  round-off  errors 
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Fig.  5   Flow  Chart  of  C0PLANAR1 


Note:  At  least  one  P  must  be  non-zero 
or  STADET  will  not  operate. 
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can  accumulate  o   For  this  reason  the  "zero1'  which  must  be 
used  in  the  program  is  not  the  mathematical  zero  but  an 
allowable  tolerance.   The  first  portion  of  this  section 
examines  the  coordinates  of  the  points  of  application  and 
determines  the  value  of  "zero"  to  be  used  throughout 
"STADET".   Appendix  L  discusses  the  method. 

Vectors  are  generated  between  the  first  point  of 
application  and  the  other  points.   These  vectors  are  ex- 
amined to  find  the  first  non-zero  one.   If  none  of  them 
is  non-zero,  there  is  only  one  point  of  support,  the  struc- 
ture is  unstable,  and  a  return  is  made  to  the  Input-Output 
Program. 

If  a  non-zero  vector  is  found,  the  subroutine  deter- 
mines whether  it  is  the  last  one  or  not.   If  it  is  not 
the  last,  the  other  vectors  are  examined  to  find  whether 
any  others  are  non-zero.   If  the  first  non-zero  vector  is 
the  last  of  the  array  of  position  vectors,  or  if  all  after 
the  first  are  zero,  the  points  of  application  of  the  links 
lie  on  one  line. 

Having  determined  that  the  points  of  application  are 
on  a  line,  the  subroutine  considers  the  number  of  links. 
A  three-dimensional  linkage  system  is  unstable  because  the 
links  intersect  a  line  and  only  a  planar  problem  has  a 
unique  solution.   Therefore,  if  there  are  more  than  three 
links,  the  system  is  redundant.   If  there  are  less  than 
three  links,  the  system  is  unstable.   If  the  number  of 
links  is  not  satisfactory,  the  appropriate  flag  is  set. 
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Where  a  second  non-zero  vector,  P£,  is  found,  its 
cross  product  with  the  first,  P-p  is  formed.   The  cross 
product  vector,  X^,  is  examined  to  find  out  if  it  is  zero. 
If  it  is  zero,  P^  and  P2  are  parallel.   Since  P^  and  P2 
have  the  same  initial  point,  this  mades  them  coincident 
and  indicates  that  the  first  three  points  of  application 
are  on  a  line.   The  subroutine  then  checks  to  see  if  all 
of  the  position  vectors  have  been  examined.   If  they  have 
not  all  been  examined,  the  search  for  another  non-zero 
vector  continues.   If  they  all  have  been  examined  and  all 
of  the  points  of  application  are  found  to  lie  on  one  line, 
an  investigation  of  the  number  of  supporting  links  is  made, 
as  above o 

If  a  cross  product  vector,  as  found  above,  is  not 
zero,  a  check  is  made  of  the  number  of  supporting  links. 
If  there  are  only  three  links,  there  are  only  two  position 
vectors  which  (since  they  are  not  parallel)  determine  a 
plane.   If  there  are  more  than  three  links,  there  are  more 
than  the  two  vectors  that  were  crossed  to  farm  X-  . ~  To 
check  whether  the  points  of  application  are  in  a  plane, 
it  is  necessary  to  dot  these  other  vectors  with  X-.  .   The 

dot  product,  X,  ,  a  scalar,  is  compared  with  zero.   If  all 
the  dot  products  are  zero,  the  points  of  application  lie 
in  one  plane. 

With  the  points  of  application  of  the  links  in  a  plane 
a  stable  system  is  possible  with  either  a  two-  or  three- 
dimensional  linkage  system.   If  all  the  links  lie  in  the 
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plane  that  contains  the  points  of  application  of  the  links , 
the  problem  is  one  in  two  dimensions.   To  check  for  a  two- 
dimensional  problem,  the  vectors  made  up  of  the  direction 
cosines  of  the  links  are  dotted  with  X-  and  the  dot  products 
are  compared  with  zero.   Since  the  link  cosine  vectors  have 
one  point  in  the  point-of-application  plane,  a  zero  dot 
product  indicates  that  the  vector  lies  entirely  in  that 
plane.   If  a  planar  problem  is  what  is  involved,  i.e.,  all 
the  dot  products  are  zero,  the  proper  number  of  links  is 
checked  for,  as  before. 

When  the  points  of  application  are  not  in  a  plane  or 
when  the  links  are  not  in  the  point-of-application  plane, 
a  three-dimensional  problem  is  involved  and  six  links  are 
required  for  stability.   Suitable  error  flags  are  set  if 
there  are  not  six  supports. 

At  this  point  it  has  been  determined  whether  a  two- 
or  three-dimensional  problem  is  involved  and  that  the 
correct  number  of  supporting  links  is  involved  for  the 
type  of  problem .   Next,  instability  due  to  incorrect 
orientation  of  the  links  is  investigated, 

b  o   PARALLELl 

This  section  checks  to  see  if  the  supporting 
links  are  parallel  or  lie  in  parallel  planes.   The  first 
link  cosine  vector  is  crossed  with  each  of  the  succeeding 
ones  and  the  resulting  cross  product  vectors  are  compared 
with  zero.   If  all  of  the  cross  products  are  zero,  all  of 
the  links  are  parallel. 
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Fig.    6      Flow    Chart    for  PAkALLELl 
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The  second  portion  of  this  section  is  concerned  with 
discovering  whether  a  situation  similar  to  the  one  shown 
in  Fig.  3  exists  or  note   In  the  above  test  a  non-zero 
cross  product  will  result  if  two  of   the  links  are  not 
parallel o   (If  the  situation  is  like  that  of  Fig,  3,  the 
cross  product  vector  is  in  the  direction  of  motion  of  the 
rigid  body0)   If  all  of  the  links  are  perpendicular  to 
this  cross  product  vector,  the  case  is  indeed  similar  to 
that  of  Figo  3»   This  section  checks  whether  the  links 
are  perpendicular  by  dotting  each  of  the  link  cosine 
vectors  with  the  cross  product  vectors  and  comparing  the 
resultant  scalar  with  zero.   If  all  of  the  dot  products 
are  zero,  all  of  the  link  cosine  vectors  are  in  parallel 
planes .   If  all  are  in  parallel  planes,  the  appropriate 
flag  is  set;   if  not,  the  program  goes  on  to  the  next 
section. 

c .  PARALLEL2 

Even  though  all  of  the  links  are  not  parallel, 
the  supporting  linkage  system  is  unstable,  by  criterion 
3.,  if  four  of  the  six  links  are  parallel.   Therefore, 
a  separate  check  is  made  for  this  possibility.   The  cross 
product  test  is  used  along  with  a  counting  system.   If  the 
linkage  system  passes  this  test,  the  remaining  check  for  a 
three-dimensional  system  is  done  by  CONCURRENTl. 

d .  PLANAR1 

This  section  of  the  subroutine  checks  the  three- 
link  cases  to  see  if  the  links  are  in  a  plane  and  if  the 
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Fig,  7  Flow  Chart  for  PARALLEL2 
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Figo  S  Flow  Chart  of  PLANAR1 
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links  are  parallel.   In  Fig*  2  there  are  two  paths  from 
C0PLANAR1  to  PLANARlo   The  first  of  these  three-link  cases 
has  the  points  of  application  on  a  line.   In  the  other 
case  the  points  of  application  are  in  a  plane  and  it  is 
known  that  the  links  are  in  the  same  plane.   For  this 
second  case  there  is  some  duplication  of  testing,  but 
the  extra  operating  time  is  negligible  and  there  was  some 
convenience  in  programming  this  way. 

This  section  crosses  the  first  link  cosine  vector 
with  the  second c   If  the  cross  product  is  zero,  the  first 
is  crossed  with  the  third  link  cosine  vector.   If  this 
cross  product  is  also  zero,  the  three  links  are  parallel 
and  the  structure  is  unstable.   If  the  second  cross  product 
is  not  zero,  it  is  dotted  with  the  second  link  cosine 
vector o   If  the  dot  product  is  zero,  all  three  link  cosine 
vectors  are  in  the  same  plane,  a  stable  planar  problem  is 
possible  and  the  program  goes  to  CONCURRENTl.   If  the  dot 
product  is  not  zero,  the  three  links  are  not  coplanar  and 
the  system  is  unstable. 

When  the  first  cross  product  is  not  zero,  it  is  dotted 
with  the  third  link  cosine  vector  and  the  dot  product  is 
compared  with  zero.   If  the  dot  product  is  not  zero,  the 
system  is  unstable,  as  above-   If  the  dot  product  is  zero, 
the  program  goes  to  C0NCURRENT1. 

e.   C0NCURRENT1 

The  remaining  check  is  to  find  if  the  lines  of 
action  of  the  links  intersect  at  a  point.   In  order  to 
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Fig.  9  Flow  Chart  of  CONCURRENT! 
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be  able  to  operate  mathematically,  it  is  necessary  to 
convert  the  link  cosine  vectors  to  lines  in  space.   The 
line  form  which  was  chosen  was  that  of  the  intersection 
of  two  planes,   (See  Appendix  G  for  discussion  of  sub- 
routine PLANE. ) 

All  of  the  planes  are  determined  and  the  coefficients 
are  stored  in  the  array  PLA .   The  equations  in  PLA  are 
examined  four  at  a  time  but  with  an  overlap,  i.e.,  equations 
one  through  four  are  examined  first ,  then  equations  three 
through  six;  and,  if  a  six-link  problem  which  gives  12 
equations  is  involved,  the  process  continues  until  all 
12  have  been  examined. 

The  three-link  system  must  have  at  least  one  inter- 
section or  the  links  would  all  be  parallel,  a  possibility 
which  has  previously  been  eliminated.   In  a  six-link 
problem  at  least  four  of  the  lines  of  action  of  the  links 
must  intersect  to  make  the  linkage  system  unstable.   There- 
fore, the  system  of  examining  pairs  of  lines  whose  equations 
are  adjacent  in  the  array  PLA  must  find  the  invalidating 
intersection  and  a  more  complicated  checking  system  is  not 
required . 

The  coefficients  of  four  equations  are  treated  as 
though  they  were  the  coefficients  of  four  simultaneous 
equations  in  four  unknowns.   They  undergo  the  matrix 
singularity  test  of  REAC1  (See  Appendix  F)„   If  the  de- 
terminant of  the  matrix  is  not  zero,  the  set  of  equations 
is  consistent  and  the  first  two  lines  do  not  intersect. 
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If  the  next  overlapping  set  of  four  is  also  consistent, 

three  lines  do  not  intersect  and  the  link  support  system 

is  satisfactory.   This  would  complete  the  checks  by  "STADET". 

On  the  other  hand,  if  the  determinant  is  zero,  the  set 
of  equations  is  dependent  or  inconsistent  and  the  two  lines 
intersect  or  are  parallel.   If  the  lines  intersect,  the 
point  of  intersection  can  be  found  by  solving  three  of  the 
equations  simultaneously.   If  the  lines  are  parallel,  the 
four  combinations  of  these  four  equations,  taken  three  at 
a  time,  are  all  inconsistent.   No  account  is  taken  of  this 
parallelism  since  it  is  already  known  that  less  than  four 
links  are  parallel.   The  next  overlapping  set  of  four  equa- 
tions is  then  considered. 

When  an  intersection  is  found,  the  coordinates  are 
substituted  into  the  equations  of  all  the  planes.   The 
number  of  equations  which  are  satisfied  is  counted.   In 
the  case  of  a  planar  problem,  if  six  are  satisfied,  all 
the  lines  of  action  intersect  and  the  linkage  system  is 
unstable.   In  the  three-dimensional  case,  if  eight  or 
more  are  satisfied,  the  linkage  system  is  unstable.  Where 
instability  exists,  the  appropriate  flag  is  set. 

This  completes  the  checks  made  by  "STADET" .   The 
rigid  body  is  supported  by  a  system  of  links  whose  re- 
actions can  be  determined  by  simple  statics  or  the  appro- 
priate flag  has  been  set  in  NOSOL. 

The  one  remaining  possible  error,  that  of  not  having 
the  load  in  the  correct  plane  for  a  planar  problem,  is  not 
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strictly  of  concern  in  deciding  whether  the  system  of 
links  is  stable  and  determinate-   The  check  for  this 
error  is  discussed  in  Appendix  C,  Input-Output  Programs 
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3c   Testing  of  Program,, 

In  order  to  devise  a  logical  system  of  tests  the  sim- 
plified flow  chart  of  Fig.  2  was  used.   This  shows  the 
questions  that  the  various  sections  of  the  subroutine 
answer o   In  all,  there  are  25  distinct  paths  that  can  be 
followed o   Only  four  of  these  paths  lead  to  a  satisfactory 
problem,  i.ea,  one  which  can  be  solved  by  simple  statics. 
Figc  10  lists  the  types  of  problems  which  were  solved .  The 
complete  print-outs  of  the  solved  problems  are  collected 
in  Appendix  K. 

Problems  26  through  29  are  the  complements  of  the 
satisfactory  problems  of  the  first  25.   The  "bar  reactions" 
output  of  TEST1A  was  used  as  the  "allowable  bar  loads"  of 
TESTlBc   The  output  of  TEST1B,  the  "allowable  load  torsor*" 
should  be,  and  is,  the  "resultant  torsor"  of  TEST1A. 

In  problems  11  and  17,  where  there  are  more  than  six 
links ,  the  storage  process  is  such  that  the  last  link  is 
stored  on  top  of  the  first*   This  results  in  a  nonsense 
print-out  of  the  structure  data.   In  problem  11,  for  ex- 
ample, the  first  link  has  two  direction  cosines  of  one. 

In  the  use  of  TEST1A  to  solve  a  planar  problem  it 
is  necessary  to  consider  the  possibility  of  having  an 
invalid  problem  because  the  load  is  not  in  the  same  plane 
as  the  structure.   Since  this  possibility  is  caused  by  an 


vTorsor  --  A  coined  word  (combination  of  torque  and 
tensor)  used  to  indicate  a  generalized  resultant  which 
is  composed  of  a  couple  and  a  force. 
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Fig=  10.  List  of  Problems 
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error  on  the  part  of  the  one  setting  up  the  problem  and 
is  not  strictly  a  factor  in  the  stability  of  the  structure 
testing  for  this  possibility  was  made  part  of  the  Input- 
Output  Program o   The  testing  of  a  planar  problem  therefore 
requires  one  step  beyond  T'STADETn.   The  subcases  of  prob- 
lems 22  and  23  show  the  results  of  such  tests. 

It  was  not  possible  to  get  an  accurate  estimate  of 
the  time  required  to  solve  a  problem .   It  takes  roughly 
one  second  to  solve  an  unsatisfactory  problem  and  about 
twice  as  long  to  solve  a  satisfactory  one„   This  differ- 
ence is  not  in  the  computation  time  but  in  the  time  re- 
quired to  write  the  output  on  the  magnetic  tape.   In  fact, 
when  a  series  of  problems  is  being  solved,  only  a  slight 
pause  can  be  detected  between  the  end  of  writing  one 
answer  and  the  beginning  of  the  output  for  the  next.  To 
these  computing  and  writing  times  must  be  added  the  time 
to  compile  the  program0   This  is  about  75  seconds  for 
TEST1A  and  about  k5   seconds  for  TEST1B, 
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4o   Conclusion o 

When  this  project  was  undertaken ,  it  was  understood 
that  this  subroutine  would  eventually  be  used  as  part  of 
a  much  more  ambitious  program 0   The  objective  of  providing 
a  routine  which  decides  whether  the  structure  may  be 
analyzed  by  the  methods  of  simple  statics  has  been  achieved. 

For  the  purposes  of  this  routine  a  decision  was  made 
to  have  storage  for  only  six  equivalent  links ,  the  maximum 
number  compatible  with  a  solution  obtainable  from  statics. 
This  is  the  only  restriction  which  must  be  removed  in 
future ,  expanded  use.   The  arrays  KPROB  (kind  of  problem) 
and  NOSOL  (no  solution)  were  set  up  to  be  used  as  keys  to 
select  the  proper  methods  of  attack .   If  N0S0L{4)  is  set 
to  one,  a  redundant  structure  is  involved  and  more  power° 
ful  methods  must  be  used  to  effect  a  solution.   KPROB 
indicates  whether  the  simplification  of  a  planar  problem 
is  valid c 

While  it  is  probable  that  a  bona  fide  planar  problem 
would  be  set  up  in  one  of  the  coordinate  planes,  it  might 
be  advisable  to  provide  a  subroutine  to  transform  the 
coordinate  system  so  that  an  oblique  plane  problem  is 
transferred  to  a  coordinate  plane „   It  may  or  may  not  be 
desirable  to  transform  the  results  back  to  the  original 
coordinate  system „ 

The  next  step  in  the  generalized  analysis  of  structures 
by  digital  computers  could  be  along  one  of  several  paths* 
Since  this  subroutine  is  concerned  with  only  the  external 
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stability  of  the  structure,  a  subroutine  to  examine  the 
internal  stability  needs  to  be  written „   The  simplifica- 
tion of  considering  only  space  frames  might  yield  a  good 
first  approximation. 

Another  path  would  be  to  either  use  a  code  or  have 
an  examining  type  of  subroutine  to  decide  the  type  of 
problem  and  then  jump  to  one  of  a  set  of  subroutines 
designed  to  solve  a  given  class  of  problems. 

It  is  felt  that  this  is  one  small  step  in  removing 
some  of  the  drudgery  from  engineering  design., 
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APPENDIX  A 
Linkage  Diagrams  of  Common  Supports 

A  link  is  a  member  capable  of  supporting  only  an 
axial  force..   It  has  spherical  hinges  at  both  ends  so 
that  it  cannot  transmit  any  moments.,   In  the  drawings 
below  the  link  is  shown  thus:   >y  ,  where  the  dot  indi- 
cates the  hinge  at  the  point  of  application.   The  other 
end  of  the  link  is  to  be  understood  to  have  a  hinge  which 
in  turn,  is  attached  to  a  frame  of  reference. 

As  discussed  in  Section  1,  Introduction,  some  of 
these  supports  are  not  the  equivalent  of  a  complete  link. 
In  these  cases  it  is  the  user's  responsibility  to  either 
modify  the  support  or  to  make  sure  that  the  reaction  is 
in  the  allowed  direction. 

Smooth  plane  support 


Because  friction  is  neglected, 
a  linkage  system  which  only 
opposes  vertical  motion  is 
obtained.   The  link  shown  may 
only  have  a  compressive  reaction 


GROUND  PLANE 


Flexible  cord 


CORD 


GROUND   The  reaction  is  along  the  cord 
PLANE   Only  a  tensile  reaction  is 
allowable. 
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Ball  and  socket 


SOCKET 
-BODY 


GROUND 
PLANE 


A  ball  and  socket  is  the 
equivalent  of  three  links 
acting  at  a  point. 


Built-in  support 


BODY 


s> 


GROUND 
PLANE 


BODY  A  built-in  support  is  the 
equivalent  of  six  links. 
The  transparent  body  is 
shown  supported  by  one  of 
many  stable  and  determinate 
arrangements  of  six  links „ 
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APPENDIX  B 
Mathematical  Tests 

The  examination  of  the  system  of  supporting  links  is 
accomplished  by  the  application  of  a  few  simple  mathemati- 
cal tests  to  the  vectors  which  are  part  of  or  are  generated 
from  the  input  data. 

Ao   The  dot  or  inner  product  of  two  vectors  A  and  B, 

neither  of  which  is  zero,  is  zero  if  and  only  if  A  and  B 

are  perpendicular,, 

A  •  B  =  |A|  |B|  cos  a     where  a   =  included  angle 

A  •  B  =  A  Rv  +  A„B__  +  A^B,,       subscripts  indicate  x, 

J    J  y,  z  components  of  A 

and  B 

Bo   The  cross  or  outer  product  of  two  vectors  A  and 

B,  neither  of  which  is  zero,  is  zero  if  and  only  if  A  and 

B  are  parallel  or  collinear., 

A  x  1  -  |A|  |^l  sin  a  H   where  u  is  a  unit  vector  per- 
pendicular to  both  A  and 
B  and  with  them  forming 
a  right-handed  set 


A  x  B  = 


i  i   k 

AAA 
Rx   Ay  Az 

B   B   B 
x  y   z 


_i,  j_  and  k  are  unit,  orth- 
ogonal base  vectors  in  the 
x,  y,  z  directions,  re- 
spectively 


Co   The  question  of  whether  of  not  lines  intersect  is 
solved  by  examining  the  equations  of  the  planes  determining 
the  lines o   If  a  set  of  four  equations  determining  a  pair 
of  lines  is  inconsistent  or  dependent,  the  lines  are  parallel 
or  intersect o   The  "matrix  singularity  test"  of  REAC1  (see 
Appendix  F)  is  used  to  find  If  the  equations  are  dependent 
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or  inconsistent o 

If  the  lines  intersect,  the  point  of  intersection 
can  be  found  by  solving  three  of  the  four  equations  simul- 
taneously. 
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APPENDIX  C 
Input-Output  Programs 

These  control  programs  provide  for  reading  in  data, 
operating  the  required  subroutines  and  printing  results. 
In  both  cases  the  first  step  is  reading  in  the  data  on 
the  equivalent  linkages .   Then  TTSTADETn  is  called  and  the 
stability  and  determinateness  of  the  support  system  is 
examined «,   The  rest  of  the  data  for  that  problem  is  then 
read  in  regardless  of  the  results  of  T!STADET?T.   The 
intermediate  operations  of  the  two  programs  will  be  con- 
sidered separately . 

TEST1A 

If  the  structure  is  stable  and  determinate,  the  re- 
sultant of  the  forces  and  couples  acting  on  the  structure 
is  calculated .   This  establishes  the  B  vector  in  the  matrix 
equation,  Ax  =  Bo   The  A  matrix,  that  array  which,  when 
used  to  premultiply  the  link  reaction  vector,  x,  gives  B, 
is  found  next.   If  there  are  six  links,  A  is  six  by  six; 
the  six  equations  represented  are  independent  and  the  values 
of  the  link  reactions  are  found  by  subroutine  REAC1. 

In  the  case  of  a  planar  problem,  "STADET"  has  examined 
the  structure  and  found  it  satisfactory,  but  it  remains  to 
be  seen  whether  the  load  is  applied  in  the  correct  plane. 
Therefore,  checks  are  made  to  see  that  the  resultant  force 
vector  is  in  the  plane  of  the  structure  and  that  the  re- 
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sultant  couple  vector  is  perpendicular  to  that  plane. 

If  there  are  only  three  links,  A  is  six  by  three; 
only  three  independent  equations  exist  among  the  six, 
and  it  is  necessary  to  solve  various  sets  of  three  equa- 
tions until  a  consistent  set  is  found  (See  Appendix  F, 
REAC1).  When  a  consistent  set  is  found,  the  results  are 
printed  out  as  discussed  below. 

TEST1B 

If  the  structure  is  stable  and  determinate,  the  A 
matrix  is  calculated.   The  allowable  loads  in  the  various 
links  are  knowns  in  this  problem.   Therefore,  the  x  vector 
is  multiplied  by  A  giving  the  allowable  load  on  the  struc- 
ture.  In  this  case  there  is  no  question  about  the  plane 
of  the  resultant  load  as  in  TEST1A. 

Output 

In  both  programs  the  intermediate  calculations  are 
suppressed  if  the  structure  is  not  stable  and  determinate. 
The  printout  then  consists  of  the  input  structure,  a 
description  of  the  problem  type  and  the  error.   If  the 
structure  is  satisfactory,  the  input  for  the  intermediate 
calculations  and  their  results  are  printed.   Since  the 
case  of  coincident  support  points  would  be  an  unusual 
blunder,  a  special  error  printout  was  not  used.   Since 
"STADET"   is  not  really  operated  in  this  case,  the  print- 
outs are: 
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NOT  A  PLANAR  PROBLEM  and 

STRUCTURE  UNSTABLE  LINES  OF  ACTION  OF  SUPPORTS 
INTERSECT  A  LINE 

(See  problem  1. ) 

Both  programs  reset  conditions  so  as  to  solve  the 
next  problem  and  then  continue  until  all  solutions  hav< 
been  made.   (See  Appendix  H  for  problem  format . ) 
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APPENDIX  D 
Subroutine  G80ST 

The  subroutine  evaluates  the  cross  product  of  two 
vectors  expressed  in  rectangular  coordinates.  A  Fortran 
expression  is  written  for  each  term  of  the  expansion  of 
CROSS. 


CROSS  = 


A 

k 

X 

A 

y 

A 
z 

35 

By 

Bz 

(AvBz»AzBv)i   +    (kzBx~kxBz)i  +    (AxBv-AyBx)k 


*y"z   "z~y 


x^y      y"X' 


=   CR0SS(l)i   +  CR0SS(2)J1  +  CR0SS(3)k 
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APPENDIX  E 
Subroutine  RES1 

The  subroutine  finds  the  resultant  of  a  system  of 
forces  and/or  couples  by  first  adding  up  the  rectangular 
components  of  either  and  then  calculating  the  effect,  if 
any,  of  the  forces  on  the  resultant  couple . 


Fx  =   i  Fxi  etc 


^=     i  Fxi  etc" 
I  =  rxl  +  ryi  +  r2k 


F  2    is  the  x  component  of  F. 


F  =  F  i  +  F  i  +  Fk 
-    x—    y^    z— 


M  =  r  x  F 


M  •  is  the  x  component  of  M± 


XI 


r  is  the  radius  vector  to  the 
point  of  application  of  the 
force 

rvS  r,.  and  r_  are  the  rectan- 
a    y       <6 

gular  components  of  r 
F  is  the  vector  force 

Fx,  F  and  Fz  are  the  rectan- 
gular components  of  F 

M  is  the  vector  moment 


Mxi  +  Myi  +  MziS 


i 
r 


x 


y 


k 
r 


Mx,  My  and  Mz  are  the  rectan- 
gular  components  of  M 
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APPENDIX  F 
Subroutine  REAC1 

This  subroutine  is  a  modification  of  the  GAUSS2  sub- 

p 
routine  used  by  C.  B.  Bailey  in  his  program  for  solving 

simultaneous  linear  equations.  Whereas  Bailey  was  solving 

the  matrix  equation,  Ax  =  B,  with  a  maximum  of  50  vectors 

B  for  each  coefficient  matrix  A,  nowhere  in  this  program 

is  more  than  one  B  vector  involved  with  a  given  A. 

The  routine  locates  the  largest  first  column  coef- 
ficient and,  if  necessary,  exchanges  the  row  containing 
this  coefficient  with  the  first  so  the  largest  element  is 
A-,  2  •   Multiples  of  row  one  are  subtracted  from  the  other 
rows  to  make  the  column  one  coefficients  below  the  first 
zero.   Next  the  column  two  coefficients  below  the  first 
are  examined  to  find  the  largest  and  rows  are  exchanged, 
if  necessary,  as  before.   This  Gaussian  elimination  process 
is  continued  until  all  the  elements  below  the  main  diagonal 
are  zero . 

At  each  reduction  step,  the  value  of  the  diagonal 
element  is  compared  with  zero.   A  very  small  diagonal 
element  indicates  ill-conditioned  equations  and  causes 
an  error  output.   This  is  the  "matrix  singularity''  test 
used  in  C0NCURRENT1.   If  the  matrix  is  non-singular,  a 
back-  solving  process  is  used  to  find  the  unknowns. 


2C.  B.  Bailey,  F2  UTEX  LINEQN,  GO  OP  Manual,  Dec,  I960. 
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APPENDIX  G 
Subroutine  PLANE 

Given  a  vector  with  direction  cosines  1,  m  and  n 
and  one   point,  the  problem  is  to  find  a  pair  of  planes 
whose  intersection  is  parallel  to  the  vector  and  passes 
through  the  point . 

First,  the  subroutine  checks  to  see  if  the  vector 
is  "nearly"  parallel  to  a  coordinate  axis.   "Nearly"  is 
defined  as  having  a  direction  cosine  greater  than  0„#, 
i.e.,  lying  within  about  37°  of  the  coordinate  directions „ 
If  the  vector  is  "nearly"  parallel  to  one  axis,  it  is 
crossed  with  vectors  along  the  perpendicular  axes . 

If  the  vector  is  not  "nearly"  parallel,  a  check  is 
made  to  see  whether  it  lies  within  about  6°  of  a  coordinate 
plane,  i.e.,  the  direction  cosine  is  less  than  0.1.   If 
it  is  close  to  a  coordinate  plane,  one  of  the  crossing 
vectors  is  chosen  perpendicular  to  that  plane. 

If  neither  of  the  above  tests  is  met,  the  given 
vector  is  crossed  with  vectors  in  the  x  and  z  directions. 

There  are  now  two  cross  product  vectors  of  the  form: 
X-l  =  Ai  +  Bj_  +  Ck 

Since   these  vectors  are  perpendicular  to  the  planes  con- 
taining the  vectors  which  formed  the  cross  products,  the 
dot  product  of  X-,  and  an  arbitrary  vector  in  one  of  the 
planes  is  zero.   An  arbitrary  vector  is  of  the  form: 
(x-x0)i  +  (y-y0)J.  +  (z-zc)k 

60 


where  (xo,yo,z  )  is  the  known  point.   The  dot  product  has 
the  form: 

A(x-xo)  +  B(y~yo)  +  C(z-zJ  =  0 
This  may  be  rearranged  to: 

Ax  +  By  +  Cz  -  D 
where:    D  =  Ax0  +  By0  +  Cz0 

The  coefficients  A,  B,  C,  and  D  are  stored  in  the  array 
PLA. 
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APPENDIX  H 

Problem  Format 


TEST1A 

Fortran  Format  Symbol(s) 

15  NPROB 
Hollerith 

15  NBAR 

3F10.S,3F14.7  VB,B 


13,12 
6F12.6 

3F14.7 


NFO,NCO 
FORCE, CORD 

COUPL 


Input  data 

Total  number  of  problems 
Problem  title  -  up  to  SO 
characters 
Number  of  links 
Direction  cosines,  coordi- 
nates of  links 
Number  of  forces,  couples 
Force  components,  points  < 
application 
Couple  components 


Repeats  from  the  Hollerith  title  for  each  problem 
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TEST1B 


Fortran  Format  Symbol(s) 

15  NPROB 
Hollerith 

15  NBAR 

3F10eS,3F14«7  VB,B 


6F12.6 


BFORCE 


Input  Data 

Total  number  of  problems 
Problem  title  -  up  to  30 
characters 
Number  of  links 
Direction  cosines,  coordi- 
nates of  links 
Allowable  link  loads 


Repeats  from  the  Hollerith  title  for  each  problem 
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APPENDIX  I 
STADET  Flags 

KPROB  =  Kind  of  Problem 

1  Support  points  are  on  a  line 

2  Support  points  are  in  a  plane 

3  Reaction  matrix  is  singular 

4  Planar  problem  (if  4  is  1) 

Not  a  planar  problem  (if  4  is  0) 

5  Resultant  load  is  not  in  plane  of  structure 

NOSOL  =  No  solution  because: 

1  Structure  unstable  -  lines  of  action  of  supports 
are  concurrent 

2  Structure  unstable  -  lines  of  action  of  supports 
are  parallel 

3  Structure  unstable  -  insufficient  number  of 
supports 

4  Structure  indeterminate  -  excess  number  of  supports 

5  Structure  unstable  -  lines  of  action  of  supports 
intersect  a  line 

6  Structure  unstable  -  supports  are  in  parallel 
planes 
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APPENDIX  J 
STADET    List 
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APPENDIX  K 
Problem  Types 

In  this  appendix  "bars"  and  ''supports"  are  used  as 
synonyms  for  "links".   See  Section  3»  Testing  of  Program, 
for  a  discussion  of  the  problems. 
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APPENDIX  L 
Definition  of  "Zero" 

The  computer  must  convert  any  non-integer  number  to 
an  approximate  binary  form  which  will  fit  its  word  length. 
The  conversion  and  length  limitation  allow  an  accuracy  of 
about  ten  decimal  digits.  With  each  calculation  round-off 
errors  accumulate.   Whereas  the  flow  charts  ask  if  a 
quantity  is  zero,  what  must  really  be  considered  is 
whether  the  quantity  is  less  than  the  maximum  expected 
round-off  error.   If  there  are  no  conflicting  requirements 
the  maximum  round-off  error  may  be  estimated  and,  with 
perhaps  some  leeway,  defined  as  "zero".   However,  where 
different  orders  of  magnitude  may  be  involved  in  the  same 
problem  or  in  subsequent  problems  or  where  it  is  advisable 
to  have  a  large  tolerance  for  one  calculation  and  a  small 
tolerance  for  another,  the  situation  becomes  far  more 
complicated.   One  cannot  define  the  "zeros"  in  terms  of 
the  number  of  significant  figures  in  an  expression;   a 
specific  magnitude  is  required. 

An  attempt  was  made  to  tie  the  size  of  the  "zero"  to 
the  size  of  the  structure.   It  was  felt  that  barring  un- 
reasonable complications,  such  as  having  a  very  small 
structure  at  a  great  distance  from  the  origin^  the  system 
of  averaging  the  coordinates  of  the  points  of  application 
and  multiplying  this  average  by  the  computer  accuracy 
would  give  a  suitable  zero.   However,  instead  of  being 
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able  to  use  10    times  the  average  coordinate,  it  was 
necessary  to  use  a  "zero"  10,000  times  as  big  to  get 
proper  answers.   Further  investigation  revealed  that 
the  difficulty  probably  involved  the  use  of  subroutine 
REAC1  in  CONCURRENTS. 

In  subroutine  REAC1  a  small  value  of  "zero"  allows 
small  diagonal  elements  which,  during  the  back-solving 
process,  generate  large  coordinates  of  the  point  of 
intersection.   When  these  large  values  are  then  used 
in  C0NCURRENT1,  the  round-off  error  becomes  much  greater 
than  the  "zero"  defined  above.   Unsuccessful  efforts  were 
made  to  modify  the  "zero"  obtained  above  in  a  logical  way 
for  use  in  C0NCURRENT1  and  REAC1. 

A  multiplier  of  10"''  gives  proper  answers,  at  least 
for  this  set  of  problems,  but  there  is  no  other  justifi- 
cation for  its  use.   The  author,  regretfully,  did  not 
have  sufficient  time  to  pursue  this  investigation. 
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